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Abstract 

The SDO/HMI instruments provide photospheric vector magnetograms with 
a high spatial and temporal resolution. Our intention is to model the coronal 
magnetic field above active regions with the help of a nonlinear force- free extrap- 
olation code. Our code is based on an optimization principle and has been tested 
extensively with semi-analytic and numeric equilibria and been applied before to 
vector magnetograms from Hinode and ground based observations. Recently we 
implemented a new version which takes measurement errors in photospheric vec- 
tor magnetograms into account. Photospheric field measurements are often due 
to measurement errors and finite nonmagnetic forces inconsistent as a boundary 
for a force-free field in the corona. In order to deal with these uncertainties, we 
developed two improvements: 1.) Preprocessing of the surface measurements in 
order to make them compatible with a force-free field 2.) The new code keeps a 
balance between the force-free constraint and deviation from the photospheric 
field measurements. Both methods contain free parameters, which have to be 
optimized for use with data from SDO/HMI. Within this work we describe the 
corresponding analysis method and evaluate the force-free equilibria by means of 
how well force-freeness and solenoidal conditions are fulfilled, the angle between 
magnetic field and electric current and by comparing projections of magnetic 
field lines with coronal images from SDO/AIA. We also compute the available 
free magnetic energy and discuss the potential infiuence of control parameters. 
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1. Introduction 



The Helioseismic and Magnetic Imager (HMI) on board of the Solar Dynamics 
Observatory (SDO) provides us with measurements based on which the photo- 
spheric magnetic field vector can be derived (jSchou et al.^ 2010). Within this 



work we describe how these measurements can be extrapolated into the solar 
corona under the assumption that the coronal magnetic field is force- free, which 
means that the Lorentz-force vanishes. We compare the resulting magnetic field 
models with observations of the coronal plasma from the Atmospheric Imaging 
Assembly (AIA), which is also onboard of SDO. 
The force-free field equations are given by 



(V X B) X B = 
V-B = 

subject to the boundary condition 

B = Bobs on the bottom boundary 



(1) 

(2) 



(3) 



where B is the 3D magnetic field and Bobs the measured magnetic field vec- 
tor in the photosphere. (Bineau, 1972 Amari, Boulmezaoud, and Aly, 2006) 



investigated the mathematical structure of these equations regarding existence, 



uniqueness and well-posedness. (Boulmezaoud and Amari, 2000) proofed the 



existence of force-free solutions for simple and multiple connected domains. (Aly, 



2005p proofed uniqueness of force-free fields for a special cylindrical configura- 



tion. Several methods have been developed to solve these equations numerically. 
For reviews see (Sakurai, 1989; Aly, 1989 Amari et al, 1997 Wiegelmann, 



2008) and within the last few years the different numerical codes have been 



intensively tested, evaluated and compared in (Schrijver et a/., 2006 Metcalf 



et fl/., 2008} Schrijver et a/., 2008). As result of a joint study (DeRosa et a/.. 



2009 ) it has been concluded that a successful application of nonlinear force free 



field (NLFFF) extrapolation methods require: 

i) Large model volumes at high resolution, which accommodate most of the 
magnetic connectivity within an active region and to its surroundings. 
The field of view of the isolated active region AR11158, as shown in Figs, [l] 
and[2j looks like a suitable candidate to fulfill this requirement. 

ii) Accommodation of measurement uncertainties in the transverse field compo- 
nent. 

This has been implemented in recent updates of different NLFFF extrapo- 



lation codes, see (Wheatland and Regnier, 2009 


Wiegelmann and Inhester, 


2010 


Amari and Aly, 2010 Wheatland and Leka, 2011 Tadesse et a/., 2011b 


). 



Preprocessing of the photospheric vector field for a realistic approximation of 
the upper-chromospheric, nearly force- free field 



■"^Preprocessing of inconsistent boundary data is in particular important for methods using the 
magnetic field vector directly as boundary condition. Grad Rubin methods use the normal mag- 
netic field and electric current (for one polarity) as boundary condition. The vertical current is 
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Figure 1. Left: Full disk SDO/HMI magnetogram, Right: Full disk AIAlTlA image. Both 
data sets have been obtained on Feb. 14. 2011 at 20:34 and have been aligned as described in 
the SDO data analysis guide [DeRosa &; Slater 2011, http:/ /www. I msal.com/sdouserguide. html . 
The rectangle outlines the sub-region (AR11158) in the vector magnetogram which is used for 
the force-free field modeling. 



As we will see in section 2.1 the HMI- vector magnetogram shown in Fig. [3] 
is almost force- free. Tools for implementing the measurement errors (previ- 
ous item) can also deal with the remaining small forces. For comparison we 
investigate also preprocessed data. 
iv) Force- free models should be compared with coronal observations. 

In Figs. [5j|6] and Table[l]we compare the force- free models with coronal images 
observed with SDO/AIA. 



2. Instrumentation and data set 



The HMI instrument (Schou et al.^ 2010) on SDO observes the full Sun at 
six wavelengths in the Fe I 6173 A absorption line. Filtergrams with a plate 
scale of 0.5^' pixels are collected and converted to observable quantities like 
Dopplergrams, continuum filtergrams, line-of-sight and vector magnetograms. 
For vector data, each set of filtergrams takes 135 seconds to be completed 
and the filtergrams are then averaged over a period of about 12 minutes. [See 
the HMI-hompage for details: http://hmi.stanford.edu/ To generate the vector 
magnetograms. Stokes parameters are derived from the averaged filtergrams and 
inverted with the help of a Milne-Eddington algorithm, an advanced version of 



the Very Fast Inversion of the Stokes Vector (VFISV) (Borrero et a/., 2010). 



The magnetic filling factor is set to be 1 in the inversion. With the amount of 



derived from the transverse magnetic field and these conditions are per construction well posed, 
even if the photospheric magnetic field vector is not force-free. Consequently preprocessing is 
not crucial for these methods. 
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Figure 2. Top left: ^los cut from full disk HMI magnetogram. Top right: Bz from vector 
magnetogram. To aline the vector magnetogram with the line of sight magnetogram and AIA 
we carried out a correlation analysis and the correlation between ^los cind Bz is 92%. Bottom: 
Same field of view seen in AIAlTlA. In the right picture we over-plotted contour lines of Bz 
(same color-code as in images above) 



information HMI provides about the line profile, this provides the most stable 
results. One could determine the filling factor in strong field regions, where it 
is expected to be close to unity. In weak field regions, however, it is difficult to 
resolve the filling factor as well as field strength. The 180° azimuthal ambiguity 
in the transverse field is resolved by an improved version of the minimum energy 



algorithm (Metcalf, 1994 Metcalf et al, 2006 Leka et al, 2009). HMI vector 



field uncertainties depend on field strength, disk position, and orbital velocity. 
Formal uncertainties due to the inversion are computed for each pixel as part 
of the normal processing. Conservatively, the random errors in the line-of-sight 
component are about 5 G, while the uncertainty in the transverse field is as much 
as 200 G in weak field regions and as little as 70 G where the field is strong. 
The zero point uncertainty in the longitudinal direction is < O.IG. Additional 
uncertainties arise because of the disambiguation and systematic errors that are 
not as well quantified. 
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Figure 3. SDO/HMI vector magnetogram observed on Feb. 14, 2011 at 20:34UT. The 
rectangular area marks the inner box, where Wf = = 1, see section |3] for details. 



Regions of interest (ROIs) containing strong magnetic fluxes are automat- 
ically identified ( [Turmon et al, 2010] ). Fig. [s] shows a vector magnetogram, 
containing AR 11158 observed on Feb. 14, 2011 at 20:34UT. After correcting 
for projection effects ( Gary and Hagyard, 199Q| ) the data have been mapped 
to a local Cartesian coordinate using Lambert equal area projection ( [Calabretta 



and Greisen, 2002). [A n overview about processing of HMI vector magnetograms 



can be found at http://Jsoc.stanford.edu/Jsocwiki/VectorMagneticField/,] For the 

magnetic field extrapolation, we bin the data to 720 km pixel (about 1^') and 
use a computational box of 300 x 300 x 160 grid points. 



2.1. Quality of the HMI vector magnetogram 



To serve as suitable lower boundary condition for a force-free modeling, vector 
magnetograms have to be approximately flux balanced and the net force and net 



torque have to vanish. (Wiegelmann, Inhester, and Sakurai, 2006) introduced 



three dimensionless parameters, the flux balance e^^^, net force balance CfQ^ce 
and net torque balance e^Qj-q^^: 



^flux 
^force 
^torque 



I Jg B^B^ 



\j,ByB^\ + \^g{Bl + Bl)-Bl 



BD-Bl 



52 

y 



!,y{{Bl + BD-Bl) 



JsVBxB^ - xByB^ 



Is x/^^ (B^ + Bl + Bl) 
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where the integrals in ^Q^ce ^torque correspond to the Maxwell stress 



tensor and it's first moment, respectively ( see ( [Molodensky, 1969 Molodensky, 
1974 Aly, 1989] )). For perfectly force-free consistent boundary conditions these 



three quantities are zero, while for real observed data this is hardly the case. For 
practical computations, however, it is sufficient these quantities become small, 
c-g- (^flux' ^force' ^torque ^ following table we list the values for the 

used HMI- data- set in the first row. All three quantities are well below unity, 
which gives us some confidence that the data might serve as suitable boundary 
condition for a force-free modeling. 



Data set 


^flux 


^force 


^torque 


HMI, Feb.l4 2011 


0.0034 


0.0564 


0.0535 


preprocessed HMI 


0.0037 


0.0002 


0.0009 


SFT Oct.26 1992 


0.0854 


0.6842 


0.8837 


Hinode Dec. 12 2006 


0.0167 


0.2727 


0.3387 


SOLIS Jun.07 2007 


0.0124 


0.6400 


0.6691 



To deal vector magnetogram data being inconsistent with the force-free assump- 
tion, we developed a preprocessing routine (Wiegelmann, Inhester, and Sakurai, 



[2006 ), which derives suitable boundary conditions for force- free modeling from 
the measured photospheric data. Applying this procedure to HMI reduces e^^^^^g 
and e^^orque further significantly (second row). It is notable that the values 
of about 0.05 in ^qj^^^q and ^torque f^^ original HMI- vector magnetogram 
(first row) is significantly lower as observed for vector magnetograms from other 
ground based and space-born missions (rows 3-5) like the Solar Flare Telescope 
(SFT), Hinode and SOLIS p| For detailed i nvestigations of these data-sets see 
(Wiegelmann, Inhester, and Sakurai, 2006 [Schrijver et a/., 2008 Thalmann, 



[Wiegelmann, and Raouafi, 2008), respectively. 



3. Nonlinear force-free field modeling 



We solve the force-free equations ([T]) - ([3|) by an optimiz ation principle as 



proposed by (Wheatland, Sturrock, and Roumeliotis, 2000) and extended by 
(Wiegelmann, 2004 Wiegelmann and Inhester, 2010) in the form: 



-i 



Wf- 



|(V X B) X B|' 



■Wd\V-Bfd^V 



^The values refer to other Active Regions and dates and the values are meant as some typical 
value-range for a particular instrument. It is planned to compare vector magnetograms for one 
particular active region and time observed with different instruments (SOLIS and HMI) and 
the corresponding force- free models (Thalmann et al., in preparation). Further investigations 
are necessary to clarify whether the good fulfillment of the force-free consistency criteria here 
is a property of this particular active region or if the HMI-measurements are more force-free 
generally. 
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+^ / (B - Bots) • W • (B - Bobs)d^S (4) 
Js 

where is a Lagrangian multiplier which controls the injection speed of the 
boundary conditions. Wf and Wd are weighting functions, which are 1 in the 
region of interest (inner 236 x 236 x 128 physical box) and drop to zero in a 32 pixel 
boundary layer towards the lateral and top boundaries of the full 300 x 300 x 160 
computational domain. W is a space-dependent diagonal matrix the element of 
which are inverse proportional to the estimated squared measurement error of 
the respective field component. In principle one could compute W from the 
measurement noise and errors obtained from the inversion of measured Stokes 
profiles to field components. Until these quantities become available, a reason- 
able assumption is that the magnetic field is measured in strong field regions 
more accurately than in the weak field and that the error in the photospheric 
transverse field is at least one order of magnitude higher as the line-of-sight 
component. Appropriate choices to optimize u and W for use with SDO/HMI 
magnetograms are investigated in this paper. For a detailed description of the 



current code implementation and tests we refer to (Wiegelmann, 2004) for the 



basic code and (Wiegelmann and Inhester, 2010) for a description and tests of 



slow boundary injection and the consideration of measurement errors. For the 
first time we combine the above described algorithm with a multiscale approach 
as described in ( Wiegelmann, 200'8| ). For this work we apply our code with a 3 



level multiscale approach to a SDO/HMI data set with 300 x 300 points in x 
and y and extrapolate 160 pixel in hight z. 

3.1. Quality of the reconstructed 3D fields 

To evaluate how well the force-free and divergence-free condition are satisfied by 
the reconstructed 3D fields, we monitor a number of expressions, which are: 

L2 = f \V-B\^d^V (6) 



'V 



(7) 



Lioo = sup |j X B| (8) 
L200 = sup|V-B| (9) 

where Li and L2 correspond to the first and second term in equation (|4]), 
respectively, with the difference that the integral is carried out in the inner 
236 X 236 X 128 physical box, where Wf = Wd = excluding the buffer boundary 
of 32 pixel towards the lateral and top boundary of the computational box. For 
potential and linear force-free fields the values correspond to the discretisation 
error. Also investigated in the inner box is the sine of the current weighted 
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average angle (Tj between the magnetic field and electric current (see ( Wheatland, 



Sturrock, and Roumeliotis, 2000 Schrijver et al.^ 2006) for details) and Lioo and 



norms for the Lorentz force and divergence, respectively. 
3.2. Code setup and choice of free parameters 

Before we perform nonlinear force-free extrapolations we use the vertical com- 
ponent Bz of the HMI-data to compute a potential and a linear force-free field 



{aL = 2.5, a = 1.16-10 m ) with a Fourier transform method ( Alissandrakis, 
[T98T| Here, a is the linear force-free parameter (V x B = aB), which is 



for a potential field. For a linear force-free field we calculate, as suggested by 



(Hagino and Sakurai, 2004), an averaged value a = jj.oJzSign{Bz)/ "^{Bzl, 
where Jz = ^ ~ vertical current in the photosphere. 

For nonlinear force-free fields we minimize the functional |4j we vary the 
Langrangian multiplier u and the Mask W, which we want to optimize. For 
cases A-H we choose W = 5T/max(5T). This seems to be a reasonable choice 
as the measurement error in the transverse field is higher in weak field regions. 
We vary the Langrangian multiplier u between 0.1 and 0.0001 for cases A-F. The 
lower the value of the slower the observed boundary becomes injected, so that 
the code has more time to relax towards a force- free state. The computing time 
increases with a power law when u decreases (Time[h] ~ O.OlSh • see 
Fig. |4|i and column ten in table [T]). The relation between u and asin(crj), can be 
be approximated also by a power law (asin(crj) [degree] 45.7° • i^^-^^), see Fig. 
^p)' The values for force- and divergence- freeness Li and L2 are slightly higher 
for u = 10~^ than for 10~^, but the general trend is that Li and L2 decrease 
with decreasing u in form of a power law, L ~ 98 • z^^-^^, with L = Li + L2, 
see Fig. It seems that the choice u = 0.001 is optimal, as higher values of 
u correspond to worse fulfilment of all force- free consistency criteria (cTj, ^1,^2) 
and a lower u only increases the computing time drastically, but does not or 
hardly improve the solution. 

In cases G and H we investigate the influence of preprocessing on the re- 
sult. We used a standard-preprocessing parameter set /ii = /i2 = 1,^3 = 
0.001, /i4 = 0.01. These parameters control the amount of force-freeness, torque- 
freeness, nearness to the actually observed data and smoothing, respectively 



(see (Wiegelmann, Inhester, and Sakurai, 2006) for details on preprocessing). 



NLFFF-extrapolations have been carried out here for u = 0.01 and v = 0.001 
and we find that Li and L2 are smaller for extrapolations from preprocessed 
data. Similar as in the unprocessed case I/i,I/2,crj decrease with decreasing z^, 
while the computing time increases. The computing time for preprocessed data 
is about a factor of a three = 0.01) or two {v = 0.001) higher compared with 
the unprocessed cases. The angle between magnetic field and current asin(crj) 
does not improve, however, and becomes worse (factor 1.3) for u = 0.001. A 
reason for the L- values becoming lower, while aj remains the same, might be 
that some strong current-peaks are smoothed out by preprocessing. If we consider 
that without preprocessing (case E) Li and L2 are already of the order of the 
discretisation error of the potential field, the lower value dj and the shorter 
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computing times, we conclude that preprocessing is not necessary for this data 
set. 

In the cases I-L we investigate the effect of different mask- functions. We 
choose a unique mask in cases I and J, which means that we do not consider 
different errors in high and low field strength regions in the photospheric vector 
magnetogram. As one can see (in comparison with cases B and E) all three 
force-free consistency criteria are worse and consequently one should not use 
a unique mask. Interesting are the cases K and L, where we choose the mask 
W = {Bt / vcidix{BT))'^ ' This choice gives more weight to strong than to weak 
regions, similar as in the linear cases A-H, but prefers strong regions significantly 
more. For u = 0.01 (case K compared with B) Li and L2 are better and cij worse. 
For u = 0.001 (case L compared with E) all three criteria are fulfilled somewhat 
better for the quadratic mask function. The computing time is, however, a factor 
of 1.6 longer for the quadratic case (L). We conclude, that if u is sufficiently low, 
than the final equilibrium is quite robust regarding the exact choice of the mask 
function profile. 

For comparison we also extrapolated the magnetogram with the old code 
version, which uses a fixed lower boundary and does not contain a mask and 
Lagrangian multiplier (cases M and N without and with preprocessing, re- 
spectively). Without preprocessing, all three force- free consistency criteria are 
fulfilled worse as for the new code, because the old code has no possibility to 
correct for inconsistencies in the magnetograms. Preprocessing improves the 
result for Li and L2 by about a factor of four, but aj hardly improves. The 
computing time for the old code (case N) is, however, significantly lower (factor 
of 5, 9 compared with E,L, respectively) as for the best runs with the new code. If 
we compare results of the old code (case N, after preprocessing) with results from 
the new code with similar computing times (cases C and G) the performance 
are similar. Higher computing times are the price we have to pay to get better 
force-free consistent equilibria. 

In column nine in Table [l] we present the ratio of the total magnetic energy to 
the energy of a potential field E/Eq. While the correct value is a priori unknown, 
this criteria cannot serve directly as a quality measure of the reconstructed 
NLFFF-fields, except that E/Eq should be greater than unity. E/Eq is an 
important quantity, as it defines an upper limit for the free magnetic energy, 
which could become converted in kinetic and thermal energy during eruptions. 
Taking the average of ah 14 NLFFF-models we find E/Eq = 1.20 ± 0.06, and if 
we consider only the best cases, say where asin(<jj) < 10° (cases D,E,F,H,J,L) 
one finds E/Eq = 1.24 ±0.03. For long time series it will probably not be possible 
to extrapolate all magnetograms with several different parameter sets, but we 
propose to do this for a some magnetograms within a time series to derive 
an error-estimation for E/Eq. We find that this quantity is not significantly 
influenced by the chosen parameter set (value of Lagrangian multiplier mask 
function profile VK, preprocessing) if the force- free consistency criteria 1/1,^2, 
are fulfilled and one should check them for each extrapolation from vector mag- 
netograms. We find that the Loo norms for Lorentz force and divergence behave 
very similar as the integral forms. Therefore both norms can alternatively be 
used to evaluate the quality of the extrapolated field. 
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Figure 4. Relation between the Lagrangian multiplier u and computing time (panel a), 
sin"-*^ aj (panel b), L = Li + L2 (panel c). Panel d shows the relation of computing time 
in hours to sin~ aj. The values correspond to NLFFF-cases A-F in Table [1] The dashed lines 
show power law fits in all panels (in panel c the lowest value u = 0.0001 has been excluded for 
the power law fit.) 



4. Comparison with AIA-images 

As pointed out in ("DeRosa et al.^ 2009") force-free field models should be com- 
pared with coronal observations in order to quantify to which extend they 
correctly reproduce the coronal magnetic field configuration. In Fig. |5]we show 
some arbitrary chosen force-free field lines (from NLFFF-model-E) in compar- 
ison with AIA-images in different wavelengths. See (Lemen et aL, 2011) for an 



overview on AIA. Qualitatively it looks like the field lines reasonably agree with 
the observed plasma-loops. Deviations, however, are also clearly recognizable. 
In the following we aim to estimate the difference between force-free field lines 
and plasma loops quantitatively. While stereoscopic reconstructed loops in 3D as 



used in (DeRosa et a/., 2009 Conlon and Gallagher, 2010) are not available, we 



compare our results with an AIA- image from one viewpoint. A basic assumption 
is that the plasma is frozen into the magnetic field, hence the plasma loops outline 
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Figure 5. AIA-images of AR11158 in different wavelengths, observed on Feb. 14, 2011 at 
20:34UT. Over-plotted are some selected field lines from NLFFF model E. 



the magnetic field lines. Consequently, the gradient of the intensity parallel to 
the magnetic field lines should be small. We are interested to detect bright loops 
(high intensity) and to quantify the deviation of projected field lines (computed 
from the NLFFF-models) and plasma loops (as visible in AIA-images) as 



C 



Jgl{s)ds 



(10) 
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Potential Field 







Figure 6. AIA-171 A images with field hnes as computed from selected models. The red, 
white and green loops here correspond to AIA-red, AIA- white and AIA-green in table ^ 



where I{s) is the intensity along a projected magnetic field line[^ This criterion 
finds loops of high I{s) with small intensity gradients along the loops (low 
V/(s)). For other possibilities to define the penalty function C see ( |Wiegelmann| 
et g/., 2QQ5| [Conlon and Gallagher, 2Q1Q[). 



For a quantitative comparison, we compute a number of field lines, origi- 
nating ±5 pixels around previously chosen locations. These locations are Pi = 
(80,30),P2 = (85,100) and P3 = (225,185). Field lines not closing within the 
extrapolation domain are excluded from the quantitative comparison and the 
ten field lines owing the lowest values of C are considered for further analysis. 
The average and standard deviation of C for these ten field lines in each region 
are displayed in the last three columns of table [l] respectively. The field lines 



"^For convenience the values in last three columns in table 
rounded to get 3-figure integer numbers. 







have been multiplied by 10^ and 
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origination from ±5 pixels around Pi,P2 and P3 are shown as red, white and 
green curves, respectively for some of the models in Fig. |6] The field lines corre- 
sponding to the potential field model, the NLFFF-models E and L (which best 
performed for the force-free consistency criteria, see last section) and a the old 
(fixed boundary, model M) NLFFF-code are shown in Fig. [gJ 

The white loops in Fig. [6] seem to agree reasonably well with AIA for all 
models, even those of the potential field model. The average penalty function 
for these loops is in the range 200-300 for all models. The red S-shaped coronal 
loops on the other hand cannot be identified with a potential or linear-force 
free model the (the linear model performing even worse than the potential one). 
All NLFFF-models show a much better agreement, with the penalty function 
for NLFFF-model (about 165-265) being in most cases less than half as large 
as for the potential field. The penalty function for the green loops is higher 
for all models, but the penalty function for the best NLFFF-models is about a 
factor of two better than for the potential field. It seems that these green loops 
are the most challenging loops to be reconstructed and the old NLFFF-code, 
which used a fixed boundary for the transverse potential field, performs only 
slightly better as the potential field. The best NLFFF-models (in the sense of 
most force-free, in particular cases E and L shown in Figure [6|, however, indicate 
the correct field topology also for the green loops. Using the penalty function C 
to evaluate the quality of the reconstruction clearly favors NLFFF-models over 
potential and linear FF models. This criterion is not sensible enough, however, to 
definitely favor one of the NLFFF-models. In the future one should consider to 
sophisticate the comparison of magnetic field models with coronal images, e.g., 
by applying different penalty functions, use loops- structures extracted from the 
images and doing comparisons in different AIA- wavelengths. For some earlier 
observation from SDO (for which vector magnetograms have not been released 
yet) one could also compare the results with images taken from vantage points 
with one or both of the STEREO-spacecraft or compare them directly with 
stereoscopic reconstructed 3D-loops. This cannot, however, be implemented as a 
standard diagnostic for NLFFF-models, as the angle between the two STEREO 
spacecraft and SDO becomes to large for stereoscopy. 

5. Conclusions and Outlook 

Within this work we carried out nonlinear force-free coronal field extrapola- 
tions of an isolated active region based on data from SDO/HMI. The vector 
magnetogram is almost perfectly flux balanced and the field of view was large 
enough to cover also the weak field surrounding the active region. Both condi- 
tions are necessary in order to carry out meaningful force- free computations. 
We also found, that the photospheric magnetogram satisfied well the force- 
free criteria. The net force and torque are considerably smaller than in earlier 
measurements of other active region fields with SET, Hinode and SOLIS. We 
do not know for sure, however, if this is a general property of HMI, or only 
true for this particular isolated active region. A comparison of an active region 
measurement with different instruments is planned. The data could be used 
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directly as boundary conditions for nonlinear force-free field computation and 
preprocessing the photospheric field was not necessary. The new code version 
incorporates errors of the measurement, in particular in the transverse field, 
and injects the boundary data slowly, controlled by a Lagrangian multiplier v. 
The error incorporation is controlled by using a mask-function, which is 1 for 
the most thrustworthy data and where one cannot thrust the data. Unless an 
exact error computation becomes available from inversion and ambiguity removal 
of the photospheric magnetic field vector, a reasonable assumption is that the 
field is measured more accurately in strong field regions and we carried out 
computations with the mask ex Bt and ex B^. For a sufficient small Lagrangian 
multiplier v = 0.001 we found that the resulting coronal fields are force and 
divergence free in the sense that the remaining residual forces are of the order 
of the discretisation error of potential and linear force- free fields. The weighted 
angle between magnetic field and electric current is about 5° —6°. The resulting 
field is almost identical for both masks, but computations with the oc B'^ mask 
take significantly longer (9h:13min instead of 5h:46min for the cx Bt mask. 
Injecting the boundary faster by choosing a higher Lagrangian multiplier (say 
V = 0.01) speeds up the computation (to half an hour), but the residual forces 
are higher and current and field are not well aligned. Inserting the boundary even 
slower (say u = 0.0001) leads to much longer computing times (more then 50h), 
but does not improve the solution. We conclude that the choice u = 0.001 and a 
mask oc Bt or oc B^ are the optimal choices for this data set. The computations 
have been carried out on one processor on a Linux-PC. Our code has been 
parallelized with Open-MP, but rather then processing a single magnetogram 
with a parallelized code, it is planned to process different magnetograms of a 
time series simultaneously. While the time cadence of HMI vector magnetograms 
is about 12min, NLFFF computations for one magnetogram on one processor 
take about 6h. Consequently the requirement is about 50 processors (for each 
active region) in order for the NLFFF-tools to catch up with the data stream 
from HMI. 

An important question is to which extend the optimum parameters for this 
data set can also be applied to other active regions from SDO/HMI. A key 
point is to monitor the consistency criteria of the magnetogram as well as the 
remaining residual forces, alignement of fields and currents in the reconstructed 
3D field. A comparison of the magnetic field model with AIA-images should 
also always be done. We used the 171 A channel here, because loops are well 
visible in this wavelength. The question, how coronal magnetic field models 
can be validated best by coronal observations should be, however, further be 
investigated. Magnetic field lines are 3D-structures and many field lines might 
not be filled with plasma and thus not visible in EUV-images. Also it is not 
trivial how/if one can use the different wavelength in AIA to validate coronal 
field models. 

A stumbling stone for AR-NLFFF models could be that other active regions 
are not so well isolated as AR11158 investigated here, but might be magnetically 
connected to other ARs and the quiet Sun. Such situations require full disk vector 
magnetograms and force-free computations in spherical geometry, as for example 



carried out from full-disk SOLIS-measurements in (Tadesse et a/., 2011a). Due to 
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it's very nature extrapolations from full disk magnetograms have to be calculated 
with a reduced spatial resolution or one has to accept significant longer comput- 
ing times. Global force-free coronal magnetic field models can also be used to 
specify the lateral boundaries for active region modelling for non- isolated ARs. 
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